Ehrenfest-time-dependent excitation gap in a chaotic Andreev billiard 
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A semiclassical theory is developed for the appearance of an excitation gap in a ballistic chaotic 
cavity connected by a point contact to a superconductor. Diffraction at the point contact is a singular 
perturbation in the limit h ^ 0, which opens up a gap -Egap in the excitation spectrum. The time 



scale ?i/-Eg 



^\nh (with a the Lyapunov exponent) is the Ehrenfest time, the characteristic 



time scale of quantum chaos. 
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The density of states in a normal metal is suppressed 
near the Fermi energy when it is brought into contact 
with a superconductor. The history of this proximity ef- 
fect goes back to the 1960's §. It was understood early 
on B that the energy range of the suppression is the 
inverse of the typical life time Tc of an electron or hole 
quasiparticle in the normal metal. This life time is finite 
(even at zero temperature) because an electron is con- 
verted into a hole by Andreev reflection at the interface 
with the superconductor ||^. The energy scale Ec = Ti/tc, 
known as the Thouless energy, is the product of the mean 
level spacing 5 in the normal metal and the dimension- 
less conductance of the contact to the superconductor. 
(For example, Ec = N5 for coupling via an A^-channel 
ballistic point contact.) The appearance of an excitation 
gap of the order of the Thouless energy is the essence of 
the traditional proximity effect. 

Some years ago it was realized that the prox- 

imity effect is essentially different if the normal metal 
becomes so small and clean that scattering by impuri- 
ties can be neglected. This applies to a quantum dot in 
a two-dimensional electron gas and because of the 
resemblance to a billiard (cf. Fig. 1) one speaks of an 
"Andreev billiard" [|ri|,^. Depending on the shape of 
the billiard, the classical dynamics varies between inte- 
grable and chaotic. No excitation gap is induced by the 
proximity effect in an integrable billiard [^01. An exci- 
tation gap does appear in a chaotic billiard , but its 
magnitude is only given by the Thouless energy if the 
chaos sets in sufficiently rapidly [^,^. 

The characteristic time scale of quantum chaos is the 
Ehrenfest time te = ln(L/AF), defined in terms of 
the Lyapunov exponent a (being the rate at which nearby 
trajectories diverge exponentially in time) and the rela- 
tive magnitude of the Fermi wave length Ap — 27: /k-p and 
a typical dimension L of the billiard jl^. Chaotic dy- 
namics requires a^^ <C Tc, but te could be either smaller 
or larger than Tc. In the regime te <C Tc the excitation 
gap is set as usual by the Thouless energy. Established 
techniques (random-matrix theory, non-linear cr-model) 
provide a complete description of this regime P,p"^-p^. 
The opposite regime te ^ Tc has no analog in the conven- 
tional proximity effect. Random-matrix theory is helpless 
and this regime has also shown a frustrating resilience to 



solution by means of the ballistic cr-model ^] . In particu- 
lar, no mechanism has yet been demonstrated to produce 
the hard gap at Ti/t-^ conjectured by Lodder and Nazarov 

i- 

Here we report an attack on this problem by an 
alternative approach, starting from the semiclassical 
Andreev approximation to the Bogoliubov-De Gennes 
(BdG) equation The limit te — > oo yields the Bohr- 
Sommerfeld approximation to the density of states 
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which is gapless (cf. Fig. |^). We have found that diffrac- 
tion at the contact with the superconductor is a singu- 
lar perturbation to pes that opens up a gap at the in- 
verse Ehrenfest time, and provides an intuitively appeal- 
ing mechanism for the gap phenomenon. 

We recall the basic equations. The electron and hole 
components u(r) and v(r) of the spinor wave function 
satisfy the BdG equation 
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which contains the single-particle Hamiltonian H = 
— + V{r) — E-p (with confining potential V) and the 
pair potential A(r) (vanishing in the normal metal and 
equal to Aq in the superconductor). The energy E is 
measured relative to the Fermi energy Ep = fcp , in units 
such that fi^ /2ra = 1. (In these units the mean level spac- 
ing 5 is related to the area A of the billiard hy 5 ~ Att/A.) 
We assume that the motion inside the billiard is ballistic 
{V = 0) and that the interface with the superconduc- 
tor is a ballistic point contact of width W ^ Xp (so 
that the number of channels N = 2M^/Af 3> 1 and the 
Thouless energy Ec — NS ^ S). We work in the regime 
Aq ^ hvp/W (which also implies Aq ^ Ec), to ensure 
that the excitation spectrum is independent of the prop- 
erties of the superconductor. 
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FIG. 1. Solid curve: Density of states pas of a chaotic An- 
dreev billiard (inset), which is gapless according to the semi- 
classical Bohr-Sommerfeld approximation (|l|). The dashed 
line indicates schematically the phenomenon that we seek to 
describe in this paper: The opening of a gap at the inverse 
Ehrenfest time as a result of diffraction at the contact with 
the superconductor. 

For a semiclassical description one substitutes (u, v) — 
{u,v)Ae^^ , with hS the action along a classical trajec- 
tory at the Fermi energy. The wave amplitude A is re- 
lated to the classical action by the continuity equation 
V • {A^VS) = 0, while S itself satisfies the Hamilton- 
Jacobi equation jVS'p — Ep — V (so that h\7S is the 
momentum along the trajectory). The BdG equation 
takes the form 

/ -2ikFds + SH A \ f ^ \ - f ^\ 

y A* 2ikFds - SH ) \ v ) ^ \v J ' 

(3) 

with SHu = ~A~iV^(^u). The derivative ds = 
kp^{\/S) ■ V is taken along the classical trajectory. The 
Andreev approximation consists in neglecting the term 
SH containing second derivatives of the slowly varying 
functions A,u,v. 

We consider a classical trajectory that starts as an elec- 
tron at a point q £ (0, W) along the interface with the 
superconductor, making an angle cf) £ (— 7r/2,7r/2) with 
the normal (cf. Fig. Q). The product b = gcos0 is the 
"impact parameter" . The trajectory returns to the inter- 
face after a path length £, and then it is retraced in the 
opposite direction as a hole . The coordinate s G (0, £) 
runs along one repetition of this trajectory. We count tra- 
jectories with measure dqdsiiKj) = dbd(j)^ corresponding 
to a uniform measure in phase space. Equivalently, we 
can sum over scattering channels n — 1,2, . . . N , related 
to by n « A^l sin0|. 



If we ignore the term 6H in Eq. (|3|) we recover the 
Bohr-Sommerfeld density of states l^-g]. Indeed, with- 
out SH the solution of the eigenvalue problem is 

with m = ±1, ±3, ±5 . . . running over positive and neg- 
ative odd integers. The path length ^ in a chaotic bil- 
liard varies in a quasi-random way upon varying the ini- 
tial conditions q and 0, with an exponential distribu- 
tion P{1) = £-iexp(-^/£). (The mean path length is 
t = AnkF/Ec ) The density of states 

/.oo 

p{E)^N d£P{£) ^ 6[E~E^{e)] (5) 

•^•^ m=l,3,5 

then evaluates to the pes ofEq. (0). 

The key assumption that will enable us to go beyond 
the Andreev and Bohr-Sommerfeld approximations is to 
assume that the amplitude A varies more slowly in space 
than the spinor components u and v, so that we can ap- 
proximate 6H by — (neglecting derivatives of A). The 
characteristic length scale for the spatial dependence 
of A is set by the smoothness of the confining poten- 
tial V, while the characteristic length scale for u, v is 
the contact width W. By assuming La ^ W we con- 
sider the case that diffraction occurs predominantly at 
the interface with the superconductor, rather than inside 
the billiard. Since A depends on the shape of the bil- 
liard, this is the regime in which we can hope to obtain 
a geometry- independent "universal" result. 

To investigate the effect of SH we restrict the dimen- 
sionality of the Hilbert space in two ways: Firstly, we 
neglect any mixing of the N scattering channels. (This 
is known to be a good approximation of the diffraction 
that occurs when a narrow constriction opens abruptly 
into a wide region ; it does not require smooth corners 
in the contact.) Secondly, since we are interested in exci- 
tation energies E ^ Ec, we include only the two lowest 
eigenstates to = ±1 of the zeroth-order solution (jj). [The 
contributions from higher levels are smaller by a factor 
exp{~Ec/2E).] We need to include both Ei and 
although the excitation spectrum contains only positive 
eigenvalues, because of the (virtual) transitions between 
these two levels induced by SH. With these restrictions 
we have for each scattering channel a one-dimensional 
eigenvalue problem. The effective Hamiltonian TioS is a 
2x2 matrix differential operator acting on functions of 
6. 

We write Ticff — Ti-o + Hi, where Ho corresponds to the 
Andreev approximation and Hi contains the diffractive 
effects. The zeroth-order term is diagonal. 

The relation between £ and b is determined by the differ- 
ential equation d£/db = g{b) exp{K£), which expresses the 
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exponential divergence of nearby trajectories (in terms of 
a Lyapunov exponent k = a/vp given as inverse length 
rather than inverse time). The pre-exponential g(b) is 
of order unity, changing sign at extrema of £{b). Upon 
integration one obtains 



K£{b) = -ln|K6| +0(1), 



(7) 



where we have shifted the origin of b such that 6 = 
corresponds to a local maximum ^max ^ ^ of £{b). [The 
logarithmic singularity is cut off at |k;&| < exp(— K^max)-] 
There is an exponentially large number Af{i) oc exp(K^ — 
i/i) of peaks around which Eq. (0) applies. 

To obtain the diffractive correction Tii, in the regime 
that 6H — — V^, we express the Lapacian in the local 
displacements ds and — e'^^db. Since these are ap- 
proximately orthogonal for ks < 1, we have 



The first term &l is a small correction to the zeroth or- 
der density of states. The second term 9^, in contrast, 
is a singular perturbation because it associates a kinetic 
energy with the variable b. The resulting zero-point mo- 
tion implies a non-zero ground state energy, and hence 
it is responsible for the opening of an excitation gap. 
Projecting TLi onto the space spanned by the two lowest 
eigenf unctions n = ±1 of Eq. (|j), and retaining only the 
leading order terms in 1/ k£, we find 



Hi = 



-3 



(9) 



The effective Hamiltonian can be brought into a more 
familiar form by the unitary transformation Tioff 
g-i<Ti77/4-^^^gifTi7r/4 (-^j^j^ ^. Y>aM\\ matrix), followed by 

the change of variable x — nb — Kb\n\Kb\ (in the range 
I a; I < 1). We work again to leading order in l/nl, and 
find 



-edl -i/hi\x\ 
i/\n\x\ edl 



K 



(10) 



This effective Hamiltonian has the same form as the BdG 
Hamiltonian (H), for a fictitious one-dimensional system 
having V — E-p and having a pair potential A (a;) that 
vanishes logarithmically c>c 1/ln \x\ at the origin (cf. Fig. 
H). The kinetic energy ed^ gives a finite excitation gap, 
even though e <g; 1. Let us now compute this gap. 

Since e^'^^'^/'^H^^e-''^^''/'* is a diagonal matrix, the 
spectrum of TLcS is given by the scalar eigenvalue prob- 
lem 



dx"^ 



\n\x\ 



E 



nkpK 



(11) 



The ground state energy is the excitation gap i?gap- To 
generate an asymptotic expansion of -Egap for small e, we 
first multiply both sides of Eq. (Ill]) by a factor Z"^ and 
then substitute x = X^/eZ. This results in 



dX2 
U{X) = 



iU 



^ = 



2Z 
IneZ 



1 - 



ZE 
nkpK 
2\n\X 

IneZ 



+ 0{lneZ)- 



(12) 
(13) 



We now choose Z such that Z^ 
the biharmonic equation 



In eZ and obtain 



{d'^ /dX'^ + 16ln\X\)^ = uj^, (14) 
LO = [ZE/irkpnf - AZ^''^ + ©(Z^^/^). (15) 



The ground state of Eq. ( |14D is at — 14.5. Substi- 
tuting in Eq. (Hj), and using Z"^/^ = | In e| - | In | In e| + 
0(1/ In e), we arrive at 



E, 



2TTkpK f ^ 3 In I lne| 



Ine 



2|lne 



8|lne 



+ 0(lne)-3/2^ 
(16) 



Only the leading order term is significant in view of the 
approximations made in Eq. (|o|). Restoring the original 
variables we have 



E,r 



ln(uF/aAF) 



(17) 



The Ehrenfest time Tp ~ a^^ ln(L/AF) contains the clas- 
sical length L — vp/a, which is of the order of the linear 
dimension of the billiard. 
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FIG. 2. Low-energy density of states p{E) of the effective 
Hamiltonian ( |lo|), r elated to p(tj) of the biharmonic equa- 
tion (0) by Eq. (|l|). The plot is for | lne| = 10 and has been 
smoothed with a Lorentzian. The inset shows the logarithmic 
pair potential appearing in HcS, the ground state of which is 
the excitation gap (dashed line). 
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The density of states p{u}) of the biharmonic equation 
( |r^ ) can be calculated numerically jlj]. The density of 
states p{E) near the gap is related to p{uj) by 

pi^E) = ^^^§^p[u; = cuo + 8\ \ne\{E/E,,^ - 1)], (18) 

-C'gap 

and is plotted in Fig. || for |lne| = 10. The factor 
M oc exp(7rfcFK/£^gap — Ec/Egs.p) counts the number of 
peaks in £(6) around which Tioff applies. The Bohr- 
Sommerfeld approximation corresponds to the large- 
uj asymptote p{uj) = exp(ijj/16). Since tj — ^ 1 im- 
plies E/Eg^p - 1 > 1/1 In e|, the width A^; ~ Eg^p/] lne| 
of the energy range above the gap in which the Bohr- 
Sommerfeld approximation breaks down is small com- 
pared to the gap itself. 

Because HeS has only a few levels in the range AE, the 
density of states p{E) oscillates strongly in this range. 
These levels are highly degenerate (by a factor Af) in 
our approximation. Tunneling between the levels will re- 
move the degeneracy and smooth the oscillations. (A 
small amount of smoothing has been inserted by hand in 
Fig. |.) We surmise that some oscillatory energy depen- 
dence will remain, but this is an aspect of the problem 
that needs further study. 

In conclusion, we have analyzed a mechanism for the 
"gap phenomenon" in the proximity effect of chaotic sys- 
tems. Diffraction at the contact with the superconductor 
is described by an effective Hamiltonian TioS that con- 
tains (1) a kinetic energy which vanishes in the classical 
limit and (2) a pair potential with a logarithmic pro- 
file. The resulting excitation gap ii^gap (being the ground 
state energy of Tieff) vanishes logarithmically as the ratio 
of the Fermi wavelength and a classical length scale (set 
by the Lyapunov exponent) goes to zero. The time scale 
h/Egiip is the Ehrenfest time, providing a manifestation 
of quantum chaos in the superconducting proximity ef- 
fect. 
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